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^ ' Abstract 

For practical wireless DS-CDMA systems, channel estimation is imperfect due to noise and interference. In this 
CJ . paper, the impact of channel estimation errors on multiuser detection (MUD) is analyzed under the framework of the 
replica method. System performance is obtained in the large system limit for optimal MUD, linear MUD and turbo 

> 

MUD, and is validated by numerical results for finite systems. 
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■ I. Introduction 

O 

O ■ Multiuser detection (MUD) lITTll can be used to mitigate multiple access interference (MAI) in direct- 



X 



sequence code division multiple access (DS-CDMA) systems, thereby substantially improving the system 
performance compared with the conventional matched filter (MF) reception. The maximum likelihood (ML) 
based optimal MUD, introduced in ifTSl . is exponentially complex in the number of users, thus being dif- 
ficult to implement in practical systems. Consequently, various suboptimal MUD algorithms have been 
proposed to effect a tradeoff between performance and computational cost. For example, linear processing 
can be applied, based on zero-forcing or minimum mean square error (MMSE) criteria, thus resulting in the 
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decorrelator ifTTll and the MMSE detector ||9|. For non-linear processing, a well known approach is decision 
feedback based interference cancellation (IC) [ 17|, which can be implemented in a parallel fashion (PIC) 
or successive fashion (SIC). It should be noted that the above algorithms are suitable for systems without 
channel codes. For channel coded CDMA systems, the turbo principle can be introduced to improve the 
performance iteratively using the decision feedback from channel decoders, resulting in turbo MUD ll20l . 
which can also be simplified using PIC Q- The decisions of channel decoders can also be fed back in 
the fashion of SIC, and it has been shown that SIC combined with MMSE MUD achieves the sum channel 
capacity ifTHI . 

It is difficult to obtain explicit expressions for the performance of most MUD algorithms in finite systems 
(Here, 'finite' means that the number of users and spreading gain are finite). In recent years, asymptotic 
analysis has been applied to obtain the performance of such systems in the large system limit, which means 
that the system size tends to infinity while keeping the system load a constant. The explicit expressions 
obtained from asymptotic analysis can provide more insight than simulation results and can be used as ap- 
proximations for finite systems. The theory of large random matrices llT2l lfT9ll has been applied to the 
asymptotic analysis of MMSE MUD, resulting in the Tse-Hanly equation 1 14|, which quantifies implicitly 
multiuser efficiency. However, this method is valid for only linear MUD and cannot be used for the analysis 
of non-linear algorithms. For ML optimal MUD, the performance is determined by the sum of many expo- 
nential terms, which is difficult to tackle with matrices. Recently, attention has been payed to the analogy 
between optimal MUD and free energy in statistical mechanics ifTOll . which has motivated researchers to 
apply mathematical tools developed in statistical mechanics to the analysis of MUD. In lfT3ll Q, the replica 
method, which was developed in the context of spin glasses theory, has been applied as a unified framework 
to both optimal and linear MUD, resulting in explicit asymptotic expressions for the corresponding bit error 
rates and spectral efficiencies. These results have been extended to turbo MUD in Q. It should be noted that 
the replica method is based on some assumptions which still require rigorous mathematical proof. However, 
the corresponding conclusions match simulation results and some known theoretical conclusions well. 
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In practical wireless communication systems, the transmitted signals experience fading. In the above 
MUD algorithms, the channel state information (CSI) is assumed to be known to the receiver. However, this 
is not a reasonable assumption since channel estimation is imperfect due to the existence of noise and inter- 
ference. Therefore, it is of interest to analyze the performance of MUD with imperfect channel estimates. 
For linear MUD, the impact of channel estimation error on detection has been studied in Q, ||2T1 and ^ 
using the theory of large random matrices. In this paper, we will apply the replica method to analyze the 
corresponding impact on optimal MUD, and then extend the results to linear or turbo MUD, under some 
assumptions on the channel estimation error. The results can be used to determine the number of training 
symbols needed for channel estimation. 

The remainder of this paper is organized as follows. The signal model is explained in Section II and 
the replica method is briefly introduced in Section III. Optimal MUD with imperfect channel estimation 
is discussed in Section IV and the results are extended to linear and turbo MUD in Section V. Simulation 
results and conclusions are given in Sections VI and VII, respectively. 

II. Signal Model 

A. Signal Model 

We consider a synchronous uplink DS-CDMA system, which operates over a frequency selective fading 
channel of order P (i.e, P is the delay spread in chip intervals). Let K denote the number of active users, N 
the spreading gain and P — ^ the system load. In this paper, our analysis is based on the large system limit, 
where K, N, P ^ oo while keeping ^ and ^ constant. 

We model the frequency selective fading channels as discrete finite-impulse-response (FIR) filters. For 
simplicity, we assume that the channel coefficients are real. The z-transform of the channel response of user 
k is given by 

p-i 

p=0 



where {gk{p)}p=o,...,p-i are the corresponding independent and identically distributed (i.i.d.) (with respect 
to both k andp) channel coefficients having variance ^. For simplicity, we consider only the case in which 
^ <^ 1. Thus we can ignore the intersymbol interference (ISI) and deal with only the portion uncontami- 
nated by ISI. 

The chip matched filter output at the Z-th chip period in a fixed symbol period can be written as 

1 ^ 

^(0 = ^ E ^^^^(0 + n(l), l^P,P + l, N, (2) 

where hk denotes the binary phase shift keying (BPSK) modulated channel symbol of user k with normalized 
power 1, {n(/)} is additive white Gaussian noise (AWGN), which satisfies E{\n{l)\'^} = a'^ ^ and {hk{l)} 
is the convolution of the spreading codes and channel coefficients: 

hkil) ^ Sk{l)*gkil), (3) 

where Sk{l) is the l-th chip of the original spreading codes of user k, which is i.i.d. with respect to both k and 

/ and takes values 1 and -1 equiprobably. We call the (7V+P-1) x 1 vector hfe = (/ife(O), hk{N+P-2))^ 

the equivalent spreading codes of user k^. Due to the assumption that ^ <^ 1, we can approximate 

by N for notational simplicity. Then the received signal in the fixed symbol period can be written in a vector 

form: 

r = -^Hh + n, (4) 

Vn 

where r = (r(P), r{N))^, H = (hi, hx) and b = (6i, bx)^. It is easy to show that 1, 
as P — > oo. Thus, we can ignore the performance loss incurred by the fluctuations of received power in the 
fading channels and consider only the impact of channel estimation error. 

B. Channel Estimation Error 

In practical wireless communication systems, the channel coefficients {fiffclOl unknown to the receiver, 
and the corresponding channel estimates {gk{l)} are imprecise due to the existence of noise and interference. 

^Note that cr^ is the noise variance, normaUzed to represent the inverse signal-to-noise ratio. 
^Superscript T denotes transposition and superscript H denotes conjugate transposition. 
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We assume that training symbol based channel estimation fT\ is applied to provide the channel estimates. 
On denoting the channel estimation error by = gk{l)—gk{l), {^6'fc(0} are jointly Gaussian distributed 

and mutually independent for sufficiently large numbers of training symbols Q. Therefore, it is reasonable 
to assume that {5gk{l)} is independent for different values of k and /. In this paper, we consider only the 
following two types of channel estimations. 

• ML channel estimation. It is well known that the ML estimation is asymptotically unbiased under 
some regulation conditions. Thus, we can assume that the estimation error 5gk{l) has zero expectation 
conditioned on gk{l), and is therefore correlated with gkil)- 

• MMSE channel estimation. An important property of the MMSE estimate, namely the conditional 
expectation E {gk{l)\Y}, where Y is the observation, is that the estimation error 5gk{l) is uncorrected 
with gk{l), and thus is biased. 

We assume that the receiver uses the imperfect channel estimates to construct the corresponding equivalent 
spreading codes, namely h^. Thus, the error of the i-th chip of is given by 

Shkii) = hk{i) - hk{i) 
p-i 

= ^Sk{i - l)Sgkil), (5) 

1=0 

from which it follows that the variance of 5hk{i) is given by = PYar{5gk{l)}. 

Fixing {6gk{l)} and considering {6sk{l)} as random variables, it is easy to show that 5hk{l) is asymptoti- 
cally Gaussian as P ^ oo by applying the central limit theorem to (5). Due to the assumption that ^ ^ 1, 
for any I, 6hk{l) is independent of most {5hk{rn)} since for any \l — m\ > P, 5hk{l) and 5hk{m) are 
mutually independent. Thus, it is reasonable to assume that the elements in 5hfc are Gaussian and mutu- 
ally independent, which substantially simplifies the analysis and will be validated with simulation results in 
Section VI. Similarly, we can assume that the elements of are mutually independent as well. 



III. Brief Review of Replica Method 



In this section, we give a brief introduction to the replica method, on which the asymptotic analysis in this 
paper is based. The details can be found in 0, ifTDI and ifT^ . 

On assuming P{bk = 1) = P{bk = —1), we consider the following ratio 



P(6, = l|r) E{b|6,=i} 



VN 



P{bk = -l\r) 



(6) 



where is a control parameter. Various MUD algorithms can be obtained using this ratio. In particular, we 
can obtain individually optimal (10), or maximum a posteriori probability (MAP), MUD (a^ = cr^), jointly 
optimal (JO), or ML, MUD (a^ = 0) and the MF (a^ = oo). 

The key point of the replica method is the computation of the free energy, which is given by 



J=-K{r,H) ^ K-'\ogZ{r,H) 



lim / P(r|i/)logZ(r,//)(ir, 

K^oo 



(7) 



where 



Z(r,ff)^5^P(b)exp 
{b} 



2a2 



1 



Hh 



and the overbar denotes the average over the randomness of the equivalent spreading codes. It should be 
noted that the second equation is based on the self-averaging assumption lfT3ll . 
To evaluate the free energy, we can use the replica method, by which we have 



J^k(j,H)= lim ( lim 



logS^ 



(8) 



where 



where bo is the same as the b in (4). However, it is difficult to find an exact physical meaning for {ba}^^^ „ . 
We can roughly consider ba to be the a-th estimates of the received binary symbols b. 

An assumption, which still lacks rigorous mathematical proof, is proposed in llT3l . which states that 
around = can be evaluated by directly using the expression of obtained for positive integers rir. 
With this assumption, we can regard as an integer when evaluating and {x^} as replicas of x. 

To exploit the asymptotic normality of -4= J2k=i ^kKk, a = 0, n^, we define variables {va}a=o,...,nr 



^^0 



'Ok, 



(9) 



f^a — ^^J2k=l^kbak, — 1,...,?!^. 

The cross-correlations of {va} are denoted by parameters {Qab}, where Qab — VM- With these defini- 
tions, we can obtain 



exp {K(3-'g {Q}) UK {Q} n dQab, 



(10) 



a<b 



where 



b K a=Q a<b 



and 



exp(6?{Q}) 



exp 



2al 



Vo{Q} 



X JJ^exp 



a=l 



Va{Q} 



dr + 0{K 



By applying Varadhan's large deviations theorem |,6J, S„^. converges to the following expression as K 



oo: 



lim K-HogEr,^ = sup^Q} {p-'g{Q}-I{Q}) 



K^oo 



(11) 



where I{Q} is the rate function of fixiQ}, which is based on an optimization over a set of parameters 

{Qab}a<b- 

Thus, the evaluation of the free energy J-'k{^, H) depends on the optimization of (11) over the parameters 
[Qab] and {Qab}, which is computationally prohibitive. This problem is tackled by the assumption of replica 

^5{x) is the Dirac delta function. 



symmetry; that is, Qoa = Qoa = -E, Va 7^ and Qab = Q, Qab = F,'ia < b, a ^ 0. Then the optimization 
of (1 1) is performed on the parameter set {m, q, E, F}. The optimal {m, q, E, F} are given by solving the 
following implicit expressions: 

m = /jg tanh [^/Fz + E^ Dz 
q = tanh^ (^VFz + E^ Dz 



E 



l+B{l-q) 



(12) 



_ l3-^B^{B-^+l-2m+q) 

^ - (l + B(l-g))2 

2 

where Dz = -^e~^ dz, Bn = A- and B = A. Then, the performance of MUD can be derived from the 
free energy, which is determined by m, g, i?, F. It is shown in (13] that the bit error rate of MUD is given 
by 

Pe = Q (4^) , (13) 



where Q{z) = Dt is the complementary Gaussian cumulative distribution function. Thus the multiple 
access system is equivalent to a single-user system operating over an AWGN channel with an equivalent 
signal-to-noise ratio (SNR) The parameters m and q are the first and second moments, respectively, of 
the soft output, 6fc = P{bk = 1) — P(6fc = — !)• When B = Bq (a^ = cr^), it is easy to check that m = q 
andE = Fusing (12). 

IV. Optimal MUD 

In this section, we discuss two types of receivers distinguished by whether or not the receiver considers 
the distribution of the channel estimation error. We denote the case of directly using the channel estimates 
for MUD by a prefix D, and the case of considering the distribution of the channel estimation error to 
compensate the corresponding impact by a prefix C. 

A. D-optimalMUD 

In this subsection, we discuss the D-optimal MUD, where the receiver applies the channel estimates 
directly to MUD and does not consider the distribution of the channel estimation error. When the equivalent 



spreading codes contain errors incurred by the channel estimation error, the corresponding free energy is 
given by 



J'K{r,H) = K-HogZ{r,H), 



(14) 



where H is the estimation of channel coefficients H and 



Z(r,^)^5^P(b)exp 



2(72 



r — 



N 



Hh 



We assume that the self-averaging assumption is also valid for SH = H — H, and thus (7) still holds with 
the corresponding S„ given by 



X 



JDo,...,Dnr a=0 

1 y/2nal Jr 



2al 



1 Y 



71 J. 

JJexp 

a=l 



2cr2 



1 ^ ^ Y 



AT 



We can apply the same methodology as in Section III to the evaluation of the free energy with imperfect 
channel estimation. The only difference is that we need to take into account the distribution of the channel 
estimation error. In a way similar to (9), we define 



1 ^ ^ 
k=l 



a = 1, rij. 



For ML channel estimation, 5hk is uncorrelated with h^, thus resulting in £^ | /i-jt^fe | = 1 and | /ife^fc | = 
1 + A^. Then we have 



^oVa = ^ J2k=l ^Okbak 



Va > 0, 



VaVb = J2k=i ^akhk, Va, 6 > 0. 
For MMSE channel estimation, 5hk is uncorrelated with hk, thus resulting in E{hkhk} 



(15) 



E{hl} 



1 — A^. Then we have 



- Z)fc=i bokKk, Va > 0, 

VaVb = — ^ Ef=i Kkhk, Va, 6 > 0. 



VoVa - ^ 



(16) 
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Thus, the free energy with imprecise channel estimation still depends on the same parameter set {m, q, E, F} 
as in Section III. An important observation is that the existence of {6hk} affects only the term Q{Q} in (10), 
and fixiQ} remains unchanged, which implies that the expressions for m and q are identical to those in 
(12). Hence, we can focus on only the computation of Q{Q}. By supposing that the assumption of replica 
symmetry is still valid, the asymptotically Gaussian random variables vq and Va can be constructed using 
expressions similar to those in [ 13 1. For ML channel estimation, we have 



Va= + A'l {zay/1 - q- ty/q) , a = l,...,nr, 
where u, t and [za] are mutually independent Gaussian random variables with zero mean and unit variance. 
With the same definitions of u, t and [za], for MMSE channel estimation, we have 



'I - Al {za^/l - q- t^) , a=l,...,n^. 
Substituting the above expressions into (10), we can obtain the following conclusions using some calculus 
similar to that of lfT3l . For ML channel estimation, the free energy is given by 

F(l-g) 



J^K (r, = j log (^cosh (yFz + E^^ Dz - Em 



(20) 



The corresponding E and F are given by 

^ l+S(l-g)(l+A^)' 

_ (l+Ag)/3-iB2(i?(7i+l-2m+(l+A2)g) 
(l+B(l-g)(l+A^))2 • 

For MMSE channel estimation, we can obtain 

Txi^.ii) = j log (cosh [yfz + E^^ Dz- Em - ~ 

1 (^ h^h A2Wi ^z.^ , B{B,' + l~il-Al)i2m-q)) \ 
- ^ + i + i,(i-g)(i-A^) I' (21) 
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and the corresponding E and F are given by 



l+B{l-g){l-A^' ^22) 
l3-^B^{l-Al){B-^+l-{l~Al){2m-q)) 
(l+B(l-g)(l-A^))2 • 

The corresponding output signal-to-interference-plus-noise-ratios (SINRs) of the ML and MMSE channel 
estimation are given by the following expressions, respectively. 

= (l + AD(-^ + /?(l-2m + (l + ADg))' ^''^ 

and 

(a2 + /?(! - (1 - A^)(2m-g))) 
Thus, we can summarize the impact of the channel estimation error on the D-optimal MUD as follows: 

• The factors t-Vt in (23) and 1 — in the numerator of (24) represent the impact of the error of the 

h 

desired user's equivalent spreading codes, which is equivalent to increasing the noise level. 

• The imperfect channel estimation also increases the variance of the residual MAI, which equals /3(1 — 
2m + (1 + A|)g) for ML channel estimation based systems and — (1 — A^)(2m — q)) for MMSE 
channel estimation based systems. 

• The equations that m = q and E = F are no longer valid when cr^ = cr^. Thus, there are no simple 
analytical expressions for obtaining the multiuser efficiency in a similar way to the Tse-Hanly equation 

ma. 

B. C-optimal MUD 

In this subsection, we consider the C-optimal MUD, where the distribution of the channel estimation error 
is exploited to compensate for the imperfection of channel estimation. For simplicity, we consider only the 
10 MUD (C-IO MUD). 
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1) ML Channel Estimation: When deriving the expressions of C-IO MUD, we consider a fixed chip 
period and drop the index of the chip period for simplicity. The conditional probability P i^{hk} |^fc|^ 
should be taken into account to attain the optimal detection. Thus, the a posteriori probability of the received 
signal r at this chip period, conditioned on the channel estimates and the transmitted symbols {6^}, is 
given by 



P r 



oc 



P r 



{hk} dbk}) P [{hk} 



hk 



K 

})n 



dhu 



(25) 



k=l 



where 



K 

P [{hk}\{hk]) =\{P [hk 



hi 



k=l 



and 



(hk\hi^ 



oc exp 



hi 

exp I -y 



It should be noted that the above two expressions are based on the assumption of normality and mutual 
independence of {6hk} in Section n.B. Then we have 



P[r 



{hk].{hk}) 



oc / exp 



V 



K 



Wp{hk\hk)dhk. 



(26) 



fe=i 



Let ri — r— Yl!k=2 hkbk, then the integral with respect to hi is given by 



oc exp 



exp 

/ 

V 



exp 



2A^ 



exp ( -y ) dhi 



/Nil+Al) 



2 (^n + {T^Jn) ) 



(27) 



where the factors common for different {bk} are ignored for simplicity 
Applying the same procedure for /z,2, /i/f, we obtain that 

/ 

{hk} , {bk} ) oc exp 



Pir 



\ 



+ 



l+A 



(28) 
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Thus the LR of lO MUD is given by 

Pih = l\r) 



Hh 



P{hk = -l|r) 



where cr^ = + jr^- Therefore, the channel estimation error is compensated for merely by changing the 



E{b|6,=-l}exp(-27 



ViV(l+A2) 



Hh 



(29) 



equivalent noise variance and scaling the channel estimate with a factor of — p 



l+A 



Similarly to the analysis in Section IV.A, we can define 



Vo = u. 1 - 



t- 



Then we can obtain the free energy, which is given by 



{za^J\ -q- ty/q) , a^l,...,nr 



(30) 



j log (cosh (yfz + E^^ Dz - Em - ^^^^ 

J_ (^ ^ B{l-q) \ , i?((i?o-^ + l)(l + AD-2m + g) \ 
2/3 I H ^(l + A^j^ l + ^l + B{l-q) J' 



(31) 



where B — 







-r^. The corresponding E and F are given by 



(l+A2+Bo(l+A2-<z))' 



(32) 



An interesting observation is that the equations m — q and E — F are recovered in this case. Also we can 
obtain the equivalent SINR, which is given by 



SINRml 



aUl + Al) + PAl + P{l-q) 



(33) 



The corresponding multiuser efficiency r) is given by solving the following Tse-Hanly style equation: 



(34) 

From (33), we can see that the impact of channel estimation error consists of three aspects, which are 
represented by the three terms in the denominator of the expression (33). The term a^{l + A^) embodies 
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the negative impact of the channel estimation error on the user being detected, which causes uncertainty in 
the equivalent spreading codes of this user and is equivalent to scaling the noise by a factor of (1 + A^). 
Besides implicitly affecting the parameter q in the third term, the channel estimation error of the interfering 
users also results in the term of /3A|; an intuitive explanation for this is that, since the output of 10 MUD 
can be regarded as the output of an interference canceller using the conditional mean estimates of all other 
users Q, the channel estimation error causes imperfection in the reconstruction of the signals of the other 
users and the variance of residual interference equals /?A^ when the decision feedback is free of errors. 

2) MMSE Channel Estimation: For MMSE channel estimation, the channel estimation error 5hk is un- 
correlated with the estimate hk. Thus, we have 



P (hk h) = P (Shk + hk 



oc exp 



{hk 



hi 



hk}' 



2A^ 



(35) 



Applying the same procedure as ML channel estimation, we can obtain the LR of 10 MUD, which is given 
by 



Pibk = 


l|r) 


P{bk = 


-l|r) 



VN 



(36) 



where the control parameter, or equivalent noise power, a'^ = a'^ + /3A^. Substituting B = ^ . a into 
(22), we have 



l+i?o(l-{l-Af)9) 
/3-ii?g(l-A^)(B-^-(2m-g)(l-A|)) 
(l+Bo(l-(l-A2)g))' 



(37) 



Similarly to the case of ML channel estimation, the equations m = q and = F are recovered as well. The 
equivalent output SINR is given by 



SINRmmse 



al + P{l-{l-Al)qy 



(38) 
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and the corresponding multiuser efficiency is given by solving the following equation: 




tanh^ 



( 




(39) 



The intuition behind (38) is similar to that of ML channel estimation. On comparing (34) and (39), an 
immediate conclusion is that the C-IO MUD is more susceptible to the error incurred by MMSE channel 
estimation than that incurred by ML channel estimation, when A| is identical for both estimators. 

V. Linear MUD and Turbo MUD 

We now turn to the consideration of linear and turbo multiuser detection. For simplicity, we discuss only 
ML channel estimation based systems in this section. MMSE channel estimation based systems can be 
analyzed in a similar way. 

A. Linear MUD 

The analysis of linear MUD can be incorporated into the framework of the replica method (for MMSE 
MUD, cr^ = cr^; for the decorrelator, cr^ — » 0) by merely regarding the channel symbols as Gaussian 
distributed random variables. The system performance is determined by the parameter set {m, g, p, F, G} 
and a group of saddle-point equations fT3!|. 

Particularly, when cr^ = (MMSE MUD), the parameters can be simplified to {q, E}, which satisfy 
q = and E = jii^^^^z^- The multiuser efficiency is determined by the Tse-Hanly equation llT4ll . 

1 ) D-MMSE MUD: Since the channel estimation error does not affect T{Q}, the parameters m, q and p 
are unchanged. With the same manipulation on Q{Q} as in Section IV, we can obtain the parameters E, F 
and G as follows: 



< 




(40) 
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2) C-MMSE MUD: Similarly to Section IV, the MMSE detector considering of the distribution of the 
channel estimation error is given by merely scaling H with a factor of ttW and changing to + -r^r- 

' h ' h 

Then, we have E = F,G = 0,m = q and p = 0. The corresponding multiuser efficiency is given implicitly 
by 

+ (4,) 

B. Turbo MUD 

1) Optimal turbo MUD: For optimal turbo MUD ll20ll . since the channel estimation error does not affect 
T{Q} when evaluating the free energy, the impact of channel estimation error is similar to the optimal MUD 
in Section IV, namely, the corresponding saddle-point equations remain the same as in L2J except that the 
parameters E and F are changed in the same way as in (20) and (32). 

2) MMSE filter based PIC: However, greater complications arise in the case of MMSE filter based PIC 
ll2m . where the MAI is cancelled with the decision feedback from channel decoders and the residual MAI is 
further suppressed with an MMSE filter. The corresponding MMSE filter is constructed with the estimated 
equivalent spreading codes |hfc| and the estimated power of the residual interference. In an unconditional 
MMSE filter, the power estimate is given by = | — 6^ j |, where bk is the soft decision feedback; 
and in a conditional MMSE filter, the power estimate is given by 1 — 6^,. However, this power estimate 



for user k is different from the true value 



h - h 



2 



since bk is unknown to the receiver, thus making the 
filter unmatched for the MAI. Hence, the analysis in ^ may overestimate the system performance since 
such power estimation errors are not considered there. Thus we need to take into account the corresponding 
power mismatch. For simplicity, we consider only unbiased power estimation. Note that this scenario can 
be applied to general cases where the received signal power is not perfectly estimated. 

For the MMSE filter based PIC, the powers of the residual interference are different for different users. 
Similarly to the analysis of unequal-power systems in ||4|, we can divide the users into a finite number (L) of 
equal-power groups, with power {Pi}i^-^ ^, estimated power and the corresponding proportion 
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{'^i}i=i L' obtain the results for any arbitrary user power distribution by letting L ^ oo. Confining 
our discussion to unbiased MAI power estimation, we normalize the MAI power such that J2l'=i ^i^i = 1 
and Q-iPi = 1. The equivalent noise variance is given by = Thus, the bit error rate of MUD is 



given by Q ^== J since the power of the desired user is unity. 
Similarly to the previous analysis, we define 

L 



1=1 kec, 
L 



^ l=\ k&Ci 

where Ci represents the set of users with power Pi. We can see that the uneven and mismatched power 
distribution does not affect the analysis of exp {Q{Q}), which incorporates the impact of channel estimation 
error. However, the rate function T{Q} is changed to 

X{Q} = sup QabQab - 5^ «z log Mfiy{Q}\ , (42) 

{Q} \a<b 1=1 / 

where 



Mfi^iQ} = - exp ^JPlPlEbo ^ 6^ + PiF ^ bah + ^^^'11 (^3) 

\ a=l a<b a=l J a=l 

in which {ba}a=i,...,nr Gaussian random variables. Similarly to [|4J, after some algebra, we can obtain 
the free energy, which is given by 



1=1 



1 / / / o^ . . ^ B (Br;^ + 1 - 2m + (1 + Al)q) \ 

- ^^iog(i + (i + A?.)(p-,)B)+ ^ ;^^(^_^)(,!,^./-'^' |- <44) 

Letting L ^ oo, we can obtain that 

m = E 



PPE 



1+P{F-G) 



p = E 



P[PPE^+2PF+l-PG) 

{l+P{F-G)Y 
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where the expectation is with respect to the joint distribution of P and P. 

For the unconditional MMSE filter, the expressions for m, q and p can be simplified to the following 
expressions, since P = E{P} = Af: 



1+A2(F-G) 

g = . (46) 

Al(^{AlY E^+2AlF+l-AlG^ 
^~ {l+Al{F-G)Y 

This implies the interesting conclusion that if the MMSE MUD based receiver regards the received powers 
of different users as being equal to the average received power, the multiuser efficiency will be identical to 
that of the corresponding equal-power system. It should be noted that the corresponding bit error rates are 
different although the multiuser efficiencies are the same. Thus, the analysis of the unconditional MMSE 
filter based PIC in [2| yields correct results. It should be noted that, for 10 MUD with binary channel 
symbols, this conclusion does not hold since the expressions for m, q and p are nonlinear inP. 

This conclusion can also be applied to frequency-flat fading channels. When the received power is per- 
fectly known, the multiuser efficiency of MMSE MUD is given by 

V + e14^\ = 1, (47) 



where the random variable P is the received power and the expectation is with respect to the distribution of 
P. When the receiver is unaware of the fading and uses equal-power MMSE MUD, the multiuser efficiency 
of this power-mismatched MMSE MUD is given by that of an equal-power system: 

c^n + E{P}r] 

Comparing (47) and (48) and applying the fact that, for any positive random variable x, E [ } < jf^^, 
we can see that this power mismatch incurs a loss in multiuser efficiency. 



VI. Simulation Results 
In this section, we provide simulation results to verify and illustrate the analysis of the preceding sections. 
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Figure 1 shows the bit error rates versus the variance of the channel estimation error for a D-IO MUD 
system with K = 10, N = 150, P = 50 and cr^ = 0.2. In this figure, 'independent' represents the case 
of equivalent spreading codes with mutually independent elements and 'convolution' represents the case 
in which the equivalent spreading codes are the convolutions of binary spreading codes and channel gains. 
From this figure, we can see that the assumption of independent elements in the equivalent spreading codes 
appears to be valid and the asymptotic results can predict the performance of finite systems fairly well. This 
figure also shows that D-IO MUD is more susceptible to the error of MMSE channel estimation than that of 
ML estimation. 

Figure 2 compares the bit error rates in D-IO and C-IO MUD systems with 13 — 0.5 and = 0.2. For 
ML channel estimation, the C-IO MUD achieves considerably better performance than the D-IO MUD. For 
MMSE channel estimation, the two lO MUD schemes attain almost the same performance. 

Figure 3 shows the bit error rates for MMSE MUD systems with the same configuration as in Fig. 3. 
Both the numerical simulations (for both independent and convolution models of the equivalent spreading 
codes) and asymptotic results are given for D-MMSE MUD, and match fairly well. Note that C-MMSE 
MUD achieves marginally better performance than D-MMSE MUD. 

Figure 4 shows the bit error rates of MMSE filter based PIC systems with the same configurations as 
in Fig.3. The decision feedback is from the channel decoder of convolutional code (23,33,37)8 when the 
input SINR is 3dB. In this figure, the theoretical and simulation results for the unconditional MMSE filter 
are represented with 'mismatched' and 'simulation', respectively; the results with the assumption that the 
residual interference power is known are represented by 'optimal'. We can observe that the optimal scheme, 
which assumes that the decision feedback error is known, achieves only marginally better performance. 

For Rayleigh flat-fading channels, the multiuser efficiency, obtained by numerical simulations, versus 
SNR is given in Fig. 5. In this figure, 'Equal Power' means the case of equal received power. For the 
case of Rayleigh distributed received power, the results of mismatched (regarding the received power as 
being equal) MMSE MUD and optimal (the received powers are known) MMSE MUD are represented 
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by 'Rayleigh-Mismatch' and 'Rayleigh', respectively. We can see that the numerical results verify our 
conclusion about the power-mismatched MMSE MUD in Section V.B. Also, the knowledge of received 
power provides marginal improvement in multiuser efficiency. 

In Fig. 6, we apply the results for C-MMSE MUD to obtain the optimal proportion a of training symbols, 
versus different coherence time M (measured in symbol periods) and system load /3, to maximize the spectral 

2 

efficiency given by (1 — a) log(l + rySNR), where SNR = 5dB, ry is determined by (41) and = We 
can see that the required proportion of training data increases with the system load and decreases with the 
coherence time. 

VII. Conclusions 

In this paper, we have discussed the impact of channel estimation error on various types of MUD algo- 
rithms in DS-CDMA systems by obtaining the asymptotic expressions of the system performance in terms 
of the channel estimation error variance. The analysis is unified under the framework of the replica method. 
The following conclusions are of particular interest: 

• The performance of MUD is more susceptible to MMSE channel estimation errors than ML ones. 

• The MUD schemes that consider the distribution of channel estimation errors can improve the sys- 
tem performance, considerably for ML channel estimation errors and marginally for MMSE channel 
estimation errors. 

• When the MMSE MUD treats different users as being received with equal power, it attains the same 
multiuser efficiency as the corresponding equal-power systems. 
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Fig. 1 

Bit error rate of D-IO MUD as a function of channel estimation error variance 
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Fig. 2 



Bit error rate of C-IO MUD as a function of channel estimation error variance 
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Bit error rate of MMSE MUD as a function of channel estimation error variance 
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Fig. 4 

Bit error rate of MMSE filter based PIC as a function of channel estimation error variance 
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Multiuser efhciency versus SNR for non-fading and Rayleigh fading systems 
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Fig. 6 



Optimal proportion of training symbols versus coherence times and system load 



